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Abstract. In this paper, the problem of Boolean operation on circular- 
arc polygons is considered. We present a concise and efficient method, in 
which the concept of related edges is introduced. Given two circular-arc 
polygons with M and N edges, respectively, the proposed method only 
takes 0((l + k) *log l + M + N) time for Boolean operation on circular-arc 
polygons, where k is the number of intersection points, I is the number 
of related edges. The tactics proposed in this paper are easily extended 
to deal with Boolean operation on traditional polygons. 



1 Introduction 

Boolean operation on polygons is a very basic operation, which has broad appli- 
cations in various disciplines such as computational geometry, computer graphic, 
CAD, GIS, motion planning |1I2I3I4I5] . Boolean operation on traditional poly- 
gons has been well investigated [6 7 8 B9I10I11I5I12I13I14I15] . a general track of 
development can be described as follows: from handling Boolean operation on 
simple polygons to dealing with complicate cases. For example, the algorithm 
proposed by Liang et al. [2] is limited to a rectangular clip polygon, and the 
algorithm proposed by Andreev |13j is limited to a convex clip polygon. Greiner 
et al. [7] proposed an algorithm that can deal with general polygon, i.e., con- 
cave polygon with hole and self-intersection. Rivero et al. [10] incorporated the 
concept of simplicial chains, which makes the treatment of special cases more 
simple. Peng et al. also adopted simplex theory, they optimized the algorithm 
proposed by Rivero et al. Later, Liu et al. [S] proposed some optimizations to 
Greiner et al. proposed algorithm. Recently, Martinez et al. [5] proposed a new 
algorithm that subdivides the edges at the intersection point, which provides 
a simple way of processing degeneracies. However, to the best of our knowl- 
edge, in existing literatures, there is few paper addressing Boolean operation on 
circular-arc polygons. 

To study this problem is meaningful since it has many applications in our real 
life. For instance, deploying sensors to ensure wireless coverage in sensor network 
is an important problem |16ll7j . In general, the sensing range of a single sensor, 
under the case no obstacle appears, is a circle. To contrary, when the obstacles 
appear, original sensing range of a single sensor should be cut off by the obstacles. 
In this case, circular-arc polygon appears. Further, assume we deployed a set of 



sensors in this kind of scenario, in order to verify the wireless coverage range 
of any two sensors (or all of sensors), Boolean operation on these circular-arc 
polygons are needed. As another example, assume there are a group of free 
rotating cameras, which are used for monitoring a supermarket. The viewing 
range of a single camera is a circle since it can freely rotate. Further, since there 
are various obstacles such as goods shelf, original visioning range of a single 
camera should be cut off by these obstacles. In this case, circular-arc polygon 
appears. When we want to verify if all places in the supermarket are monitored 
by these cameras, Boolean operation on these circular-arc polygons are needed. 
Last but not least, assume there are a group of free moving robots, which are 
used for guiding the visitors in a museum. Since the energy of a single robot is 
limited, the region it can move is restricted in a circle. In addition, since the 
impact of various obstacles such as exhibits, original region it can move should 
be cut off by these obstacles. In this case, circular-arc polygon appears. When 
we verify if every place in the museum can be served by at least a robot, then 
Boolean operation on these circular-arc polygons are needed. In fact, there are 
many applications in our real life just like the three examples above, for saving 
space, we do not further enumerate them. 

At hrst sight, Boolean operation on circular-arc polygons is very simple, 
we just need to calculate the intersection points, then to obtain the resultant 
polygon based on a proper traversing rule. However, as discussed in subsection 
12. 31 the efficiencies of the naive methods are unsatisfactory though they also can 
achieve Boolean operation on circular-arc polygons. 

As we know, circular-arc polygon is a special case of conic polygon, Boolean 
operation on conic polygons has been investigated though only a few papers can 
be found. For instance, Berberich et al. [TS] suggested using plane sweep method 
to achieve Boolean operation on conic polygons, but few details were discussed. 
In their paper, they focused on the arrangement of conic arcs. Recently, Gong 
et al. [IH] give detailed discussion about Boolean operation on conic polygons, 
they presented an 0(M *N) time algorithm, where M and N are the number of 
edges in the two conic polygons, respectively. No doubt, their proposed approach 
also can deal with Boolean operation on circular-arc polygons. However, it has 
some obvious drawbacks under the scenario where only circular-arc polygons and 
traditional polygons appear. One of drawbacks is the redundant and complex 
data representation, another drawback is its complicate implementation as well 
as high time complexity. 

Motivated by the reasons discussed above, we investigate the problem of 
Boolean operation on circular-arc polygons. By well considering the nature of 
this problem, we propose an concise and efficient approach, in which the classical 
plane sweep method is incorporated, more importantly, we present several new 
ideas. The novelty in our proposed method has three points: (1) we introduce 
the concept of related edges; (2) we adopt two extra sequence lists in which 
the processed related edges, decomposed arcs and intersection points are well 
organized; (3) we assign two labels to each processed related edges before it is 
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placed into red-black tree, each label has its special function. To sum up, we 
make the following contributions. 

— We devise a simple and easy to operate data structure which is convenient 
for Boolean operation on circular-arc polygons. 

— We develop a concise and efficient algorithm for Boolean operation on circular- 
arc polygons. 

— We analyse the computational complexity of our proposed algorithm, it 
shows a tighter upper bound. 

— We state that the tactics used in this paper also are effective for Boolean 
operation on traditional polygons, and the 0((l + k) * \ogl + M + N) time 
complexity also holds for Boolean operation on tradition polygons. 

The remainder of the paper is organized as follows. In Sectionj^J we introduce 
the preliminaries. We present our solution in Section [3j and draw a comparison 
with state-of-the art results in Section |4l Finally, we conclude this paper and 
present our future research directions in Section [5j 



2 Preliminaries 

In this section, we first discuss the data structure for circular-arc polygon in 
subsection 12. II Next, in subsection 12. 2[ we review standard plane sweep method 
[20 21 which will be employed in our proposed method. For ease of understand- 
ing our proposed method, we discuss two naive approaches in subsection 12. 31 In 
the rest of paper, if without stated otherwise, we use intersections and intersec- 
tion points interchangeably. 



2.1 Data structure 

Definition 1. (Circular-arc Polygon) A polygon is a circular-arc polygon such 
that the boundary of this polygon is consisting of circular arcs, or the blend of 
straight line segments and circular arcs. 

As we know, the traditional polygon can be represented by a series of vertexes. 
However, for polygon with circular arcs, this method is invalid since two vertexes 
in here cannot exactly determine a circular arc segment. Fortunately, we observe 
that two vertexes incorporated with an appendix point can exactly determine 
a circular arc segment. Taking the circular-arc polygon in Figure 1(a) as an 
example, obviously, we cannot exactly represent it if we only use five vertexes, 
{Vi, V2, V3, V4, V5}. The reason is that the circular arc V1V5 is ambiguity, it may 
be a major arc or a minor arc. Note that, when we add an appendix point V§, 
then, the ambiguity can be eliminated, namely, the circular arc ViVqVs is without 
ambiguity. 

In view of the above observation, we devise a data structure that can repre- 
sent circular-arc polygon and traditional polygon in an unified manner, thus this 
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data structure should be convenient not only for Boolean operation on circular- 
arc polygons, but also for Boolean operation between traditional polygon and 
circular-arc polygon. Specifically, we adopt a double linked list, in each node, 
there are several domains as follows. 

— Data: (x,y), the coordinate of point is stored in this domain. 

— Tag: Boolean, this domain indicates whether a point is traditional vertex or 
an appendix point. 

— Crossing: Boolean, this domain indicates whether this point is an intersec- 
tion point. 

— EE: Boolean, this domain indicates what property (entry or exit) an inter- 
section point has. 

— PPointer: Node pointer, this domain points to the previous node. 

— NPointer: Node pointer, this domain points to the next node. 

For clearness, in the rest of this paper, an appendix point is denoted by Vi, 
and a traditional vertex is denoted by Vi. For example, the circular-arc polygon 



in Figure 1(a) can be denoted by {V5, V6,Vi, V2, V3, V4}. 



2.2 Plane sweep method 

One elegant method for reporting the intersections among a set of segments is the 
plane sweep method [20121) . Their basic idea is as follows. Sorting the endpoints 
of all segments based on their x-coordinates, and put them into a heap Q. The 
plane is swept (from left to right) with a vertical line. At each endpoint during 
this sweep, the segments which intersect the sweep line are stored, ordered from 
bottom to top, in a red-black tree R. Then, to check if newly adjacent segments 
intersect each other whenever the R is modified. In this way, all intersections 
can be obtained Q. 

1 Note that, in some cases, some line segments may be vertical line segments, and many line seg- 
ments possible intersect at one point, for these problems, please refer to the literature 1201 for 
more details. 
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Taking the three line segments 11 L\, Li and L3" in Figure 1(b) as an example, 
on the second sweep line, Li is inserted into the red-black tree R, then check if 
Li intersects with its adjacent segment. In here, L\ is the only adjacent segment, 
and it intersects with Li, so the intersection point I\ is obtained. At the same 
time, we insert 1\ into Q (in x-order). On the third sweep line, L3 is inserted 
into R, the order in R is "L3, L%, Li" , no adjacent segments intersect with L3, 
we do nothing. On the fourth sweep line, it is an intersection point of L\ and 
Li. In this case, the locations of L\ and Li (in R) are swapped. After this, the 
order in R is tl L%, L\, Li" now. For L\ and L2, L\ is the upper segment, Li is 
the lower segment, we check if L\ intersects with the segment immediately above 
it, and check if Li intersects with the segment immediately below it. In here, L3 
is immediately above L\, and it intersects with L\, then, the intersection point 
Ii is obtained. At the same time, we insert Ii into Q (in x-order). Note that, 
there is no segment below Li, we do nothing. For the rest of sweep lines, similar 
processing is conducted. 

Up to now, we assume the reader have understood the idea of plane sweep 
method, if any question, please refer to paper (e.g., [20]) for more details. 



2.3 Naive approaches 

For ease of understanding our proposed method, in this subsection, we introduce 
the naive approaches. Given two circular-arc polygons CAP\ and CAPi, which 
have a number M and N of edges, respectively. For achieving Boolean operation 
on them, there are two naive approaches. 



Approach 1 For each edge in CAP\, we check if it intersects with the edges 
of CAPi. If so, we compute the intersection points, and insert them into corre- 
sponding edges. In this way, until all edges in CAPi are dealt with. As a result, 
all intersection points are inserted into the two circular-arc polygons. Next, we 
assign entry-exit properties |9I7] to these intersection points. At last, we get the 
resultant polygon by traversing the two circular-arc polygons. 

However, this method bears a heavily time cost especially when the number 
of edges in the two circular-arc polygons are much more. In this approach, since 
obtaining the intersections takes 0(M * N) time, assign entry-exit properties 
and traversing takes linear time. Thus, the total time complexity is 0(M * N). 



Approach 2 The second naive approach is using the standard plane sweep 
method. Before we discuss the specific steps of the second naive approach, we 
first give the definition of non-x-monotone circular arc and x-monotone circular 
arc. 

Definition 2. (Non-x-monotone Circular Arc) Given any circular arc, it is a 
non-x-monotone circular arc such that there exist at least one vertical line that 
intersects with the circular arc at two points. 
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Definition 3. (X-monotone Circular Arc) A x-monotone circular arc is such a 
circular arc that is at most one intersection with any vertical line. 

Obviously, if there exist non- x-monotone circulars in CAP\ or CAP2, it will 
make the standard plane sweep method invalid. So, we have to decompose them 
before executing standard plane sweep method. 

Theorem 1. (Decomposing Non-x-monotone Circular Arc) Given any non-x- 
monotone circular arc NMC , and termed its corresponding circle as C , a hor- 
izontal line I, which passes through the center of C , can decompose the NMC 
into at least two and at most three x-monotone arcs. 



Proof Since the I passes through the center of C, it must intersects the C at 
two points, say p\ and p 2 . For the NM C, it only has two cases. (1) both p\ and 
P2 are on the NMC. (2) only p\ (or p 2 ) is on the NMC. For case 1, based on 
analysis geometry, it is obvious that NMC will be decomposed into two arcs by 
I when pi (or p 2 ) is the endpoint of NMC (see Figure 2(a) I, otherwise, NMC 
will be decomposed into three arcs by I (see Figure 2(b)). For case 2, it also is 
obvious that NMC will be decomposed into two arcs by I (see Figure 2(c) I. 

In addition, since both the upper semi-circle and below semi-circle of C are x- 
monotone arcs, and any sub-arc of upper or below semi-circle also is x-monotone 
arc. Putting all together, thus the theorem holds. □ 




(a) (b) (c) 

Fig. 2. Illustration of Theorem [T] 



Steps: For each edge in the two circular-arc polygons, we check if it is non-x- 
monotone arc, if so, we decompose it into x-monotone arcs. Next, we compute all 
the intersection points based on the standard plane sweep method. Note that, 
when we use standard plane sweep method, the vertexes of CAP\ and CAP 2 
also are reported as the intersection points. So, we next have to choose real 
intersection points from all reported points. For each real intersection point, we 
insert it into the two circular-arc polygons. Once all real intersection points are 
inserted, we assign entry-exit properties |9I7) to these intersection points. At 
last, we get the resultant polygon by traversing the two circular-arc polygons. 
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Cost analysis: In this approach, decomposing a non-x-monotone arc takes 
constant time, at worst case, each edge is a non-x-monotone arc, so decomposing 
non-x-monotone arcs takes 0(M + N) time. Based on Theorem [TJ there are at 
most 3* (M+N) edges after decomposing all non-x-monotone arcs. Without loss 
of generality, we assume there are k real intersection points, so the number of all 
reported points is at most /c + 6* (M + N), thus, executing standard plane sweep 
algorithm takes 0( ( 3 * (M + N) + (k + 6 * (M + TV)) ) * log(3 * (M + N)) ) 
time [5D], i.e., 0((M + N + k) log(M + N)) time. Choosing the k real intersection 
points from all the k + Q * (M + N) reported points can be reduce to the problem 
of removing 6 * (M + N) balls from 6 * (M+N) + k different coloured balls, which 
takes 0( (6*(M+N))*(6*(M+N) + k) ) time, i.e., 0( (M + N)*(M + N + k) ) 
time. In addition, inserting these intersection points takes 0( k* (M + N) ) time, 
and traversing takes linear time. Thus, the total time is 0( (M + N + k) * (M + 
N + log(Af + AT)) )■ 

It is easy to know, in this approach, the following steps take much time: (1) 
executing the standard plane sweep algorithm; (2) choosing the real intersection 
points from all reported points; (3) inserting the real intersection points into the 
two circular-arc polygons. 

Just now, we introduce a pair of naive approaches, both of them can achieve 
Boolean operation on circular-arc polygons, but their efficiencies are unsatisfac- 
tory. In next Subsection, we present our solution. 

3 The proposed method 

As we know, if the two circular-arc polygons are disjointed, or one is totally 
included in another, the Boolean operation is very simple. For presentation sim- 
plicity, we omit the details. In the rest of the paper, we mainly focus on the case 
where the two circular-arc polygons intersect each other. 

In this section, we discuss the proposed method, which incorporates the idea 
of plane sweep method, at the same time, we propose several new ideas. On the 
whole, the proposed method can be summarized in three steps. 

1. Constructing two sequence lists in which the intersection points are well 
organized, and the consecutive decomposed arcs are well arranged. 

2. Constructing the modified polygons based on the above two sequence lists. 
In the two modified polygons, the intersection points, appendix points and 
original vertexes have been arranged, the decomposed arcs have been merged. 

3. Obtaining the resultant polygon by traversing the two modified polygons. 

3.1 constructing two sequence lists 

For ease of understanding our ideas, in this subsection, we first describe several 
important observations, Next, we describe several new tactics. At last, we present 
the detailed algorithm for constructing two sequence lists. 
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Observations There are several important observations, which motivate us to 
develop more efficient algorithm. 

— Observation 1: Vertexes from circular- arc polygon do not need to be reported 
as intersection points. 

— Observation 2: The newly adjacent segments (in red-black tree R) do not 
need to be checked (whatever they intersect or not) if they are from same 
circular-arc polygon. 

— Observation 3: In some cases, not all edges are needed in order to obtain the 
intersection points. 

— Observation 4: The obtained intersection points should be well arranged for 
facilitating the follow-up step, i.e., constructing the modified polygons. 

Observation 1 and observation 2 are intuitionistic, but the importance must 
not be ignored. Regarding to observation 3, let us take the two circular- arc 



polygons in Figure 3(a) as an example, though there are many edges in CAP\ 
and CAP2, only 4 edges (i.e., ab, be, AB, BC) are needed in order to obtain 
the intersections, we call these edges as related edges that will be defined later. 
Observation 4 is derived from the following fact. 

FACT 1. Given two circular-arc polygons with M and N edges, respectively, 
and given we have obtained the number k of intersection points, it takes 0( k * 
[M + N) ) time to insert these intersections into the two circular-arc polygons, 
if we use naive method. 



Tactics In view of those observations discussed in previous subsection, we in- 
troduce several new tactics. 

— Tactic 1: Choosing related edges before doing others. 

— Tactic 2: Constructing two sequence lists for managing the intersection points 
and the decomposed arcs. 

— Tactic 3: Assigning two labels to the edge before it is inserted into red-black 
tree. 

Next, we first introduce how to choose the related edges from two circular-arc 
polygons, followed by explaining the details of two sequence lists as well as the 
functions of two labels. 

Related edges: The true purpose of choosing related edges is to avoid oper- 
ations that are irrelevant with obtaining the final result as much as possible. 
Before presenting the details of related edges, we first give several definitions, 
which are the basis of the concept of related edges. 

Definition 4. (Extended Boundary Lines) Given a circular-arc polygon, we al- 
ways can find a y-axis (or x- axis) paralleled rectangle box, i.e., so called minimum 
bounding rectangle (MBR)Q, which just enclose the circular-arc polygon. With- 
out loss of generality, assume the coordinate of left-bottom corner of the MBR 

^ See http://en.wikipedia.org/wiki/Minimum_bounding_rectangle 
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is "(xx,yi)", the coordinate of right-top corner of the MBR is "(x2,y2)"- Then, 
the following four lines, "X—xi", "X=X2", "Y—yi", "Y =yi" are the left, right, 
bottom and top extended boundary lines of this circular-arc polygon, respectively. 

Definition 5. (Effective Axis) Given two circular- arc polygons CAP\ andC 'AP2, 
assume their MBRs are MBR\ and MBR2, the intersection set of MBR\ and 
MBR2 also is a rectangle, say ISMM . If the horizontal span of IS MM is larger 
or equal to its vertical span, then y-axis is the effective axis. Otherwise, x-axis 
is the effective axis. 

Definition 6. (Related Edges) Given two circular-arc polygons CAP\ andC 'AP2, 
we use L\, R\, T\ and B\ to denote the left, right, top and bottom extended 
boundary lines of CAP\, respectively. Similarly, we use L2, R2, T2 and B2 to 
denote the left, right, top and bottom extended boundary lines of CAP2 , respec- 
tively. 

— If the effective axis is x-axis, without loss of generality, assume that L\ < 
L2 < R\ < R2, where L\ < L2 denotes the L\ is in the left 0/L2. Then, the 
following edges are related edges. (I) edges that intersect L2 or R\; or (2) 
edges that locate between L2 and R\ . 

— To contrary, if the effective axis is y-axis, without loss of generality, assume 
that B\ < B2 < T\ < T2, where B\ < B2 denotes the B\ is in the bottom of 
B2- Then, the following edges are related edges. (1) edges that intersect B2 
or T\; or (2) edges that locate between B2 and T\. 

For example, two big rectangle denote the MBRs of CAP\ and CAP2, the 



grey rectangle denotes the intersection set of two MBRs, as shown in Figure 3(a) 
In this example, the horizontal span of grey rectangle is less than its vertical 
span, then, the x-axis is effective axis. Further, four dotted vertical lines L\, 
L2 and i?2 denote the left, right extended boundary lines of the two circular- 
arc polygons, respectively. Since edges ab and be intersect with L 2 , so they are 
related edges. Similarly, edges AB and BC intersect with R\, so they also are 
related edges. 
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Note that, in order to let the chosen related edges be well organized, when 
we choose the related edges, we follow a rule as follows. 



— If the effective axis is x-axis. 

• For the left circular-arc polygon, we detect each edge, starting from the 
leftmost edge and along the direction of counter-clockwise, to see if it is 
a related edge. 

• For the right circular-arc polygon, we detect each edge, starting from the 
rightmost edge and along the direction of counter-clockwise, to see if it 
is a related edge. 

— If the effective axis is y-axis. 

• For the top circular-arc polygon, we detect each edge, starting from the 
topmost edge and along the direction of counter-clockwise, to sec if it is 
a related edge. 

• For the bottom circular-arc polygon, we detect each edge, starting from 
the bottommost edge and along the direction of counter-clockwise, to 
see if it is a related edge. 

Definition 7. (Processed Related Edges) Given there are a number I of related 
edges, if there are non-x-monotone arcs among these related edges, we decompose 
them. Without loss of generality, assume there are a total number I + I* of edges 
after decomposing. We call them as processed related edges. 

Lemma 1. Given a number I of related edges, if there is no non-x-monotone arc 
among them, the number of processed related edges is I . Otherwise, the number 
of processed related edges is larger than I and no more than 3 * I. 

Proof When there is no non-x-monotone arc, we need not decompose any re- 
lated edge, obviously, \\related edges\\= \\processed related edges\\. 

When there exist non-x-monotone arcs, based on Theorem [TJ we have that 
\\processed related edges\\ < 3*\\related edges\\. In addition, assume there is only 
one non-x-monotone arc, based on Theorem[TJ we have that \\processed related edges j| 
> \\related edges\\ + 1. Thus, I < \\processed related edges\\ < 3 * I □ 

Two sequence lists: In the two sequence lists, each item is a compound struc- 
ture. The compound structure consists of three domains: (1) the processed re- 
lated edge; (2) the locations of intersections (if exist) on this edge; (3) a tri-value 
switch. 



iteml 
itme2 



(ab) | 1 0.7.0.9 | |~0~| 



iMll 0.1,0.5 I QD 



iteml 
itme2 



S 2 

HABT1 | 0.8,0.91 [U 



RbcTI I 0.1,0.3"| Qj] 



Fig. 4. Example of Sequence List 
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The processed related edges in each sequence list are stored along the direc- 
tion of counter-clockwise w.r.t. the original circular-arc polygon. By doing so, it 
is convenient for the follow-up step, i.e., constructing the modified polygons. As 
an example, regarding to the two circular-arc polygons in Figure [3j we build two 
sequence lists, Si and 5*2, as shown in Figure H) Edge ab is stored in Si as the 
first item, since it is prior to the edge be in the CAP\ in terms of the direction 
of counter-clockwise. Similarly, edge AB is prior to the edge BC in the CAP2 
in terms of the direction of counter-clockwise, thus, edge AB is stored in S2 as 
the first item. 

Regarding to the second domain, we use a series of decimals to denote the 
locations of intersections, e.g., a decimal 0.5 denotes that an intersection point 
just locates at the middle point of an edge. Let us revisit Figure 21 the two 
decimals 0.7 and 0.9 denote that there are two intersections on edge ab, and 
^(aiVlf = 0-?i ^a' 7 fr 2 )|f =0-9; where ||(a,7i)|| denotes the length from endpoint a 
to the intersection point I\ (cf. Figure [3(b)) , ||(a, 7 2 )|| and ||(a, &)|| have similar 
meanings. Certainly, if an edge is arc rather than line segment, then, the decimal 
denotes the length of arc (from stating point to intersection point). 

Before we discuss the functions of the tri- value switch (i.e., the third domain), 
we first give a definition as follows. 

Definition 8. (Decomposed Circular Arc) Given two circular-arc polygon CAP\ 
and CAP2. Among all the related edges, assume there is a non-x-monotone cir- 
cular arc, say NMC . We decompose it into several x-monotone circular arcs, say 
XMCi, where i € {1,2,3}. Then, each x-monotone arc XMCi is a decomposed 
circular arc. 

The value in the third domain either is "0" , "1" , or "2" . When the value is 
"0", it indicates this edge is not a decomposed circular arc. In above example, 
all related edges are line segments rather than non-x-monotone arcs, so there 
is no decomposed arc. Thus, the third domain of each item is filled with "0" , 
as shown in Figure 2) When the value is "1" or "2" , it indicates this edge is a 
decomposed arc. The reason we use two different values is to help us differentiate 
the decomposed arcs, which are from two different non-x-monotone arcs. The 
goal of differentiating these decomposed arcs is convenient for merging these 
decomposed arcs correctly and efficiently in future. 

Then, given a series of decomposed arcs, what value should be assigned to 
each single decomposed arc? The assignment rules are as follows. Specifically, 
we assign "1" to each decomposed arc, which is from the first non-x-monotone 
arc, and assign "2" to each decomposed arc, which is from the second non-x- 
monotone arc. If there are still non-x-monotone arcs, we assign "1" and "2" to 
their corresponding decomposed arcs alternately. 

For example, consider the two circular-arc polygons in Figure 5(a)[ according 



to the definition of related edges, the following edges in CAP\ are related edges: 
" ab, bgc, chd, die, ef ". Note that, since edges bgc, chd and die are non-x- 
monotone arcs, we decompose them based on Theorem[TJ e.g., bgc is decomposed 
into segments bmj and jgc, similarly, edges chd and die also are decomposed, as 
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(a) (b) 
Fig. 5. Example of Consecutive Non-x-monotone Circular Arcs 



shown in Figure [5 (b)| In this Figure, j, k and I are the decomposing points, m, 
n and o are new appendix points. After decomposing, there are a number eight 
of processed related edges that will be put into sequence list S\. So, the number 
of items in sequence list Si is eight. Specifically, segment ab as the first item, 
segment bmj as the second item, and so on. Note that, for the number eight of 
items in sequence list Si, based on the assignment rules, the values of the third 
domains should be "0", "1", "1", "2", "2", "1", "1", "0", respectively. 

Another thing we have to mention is that, we adopted sequence list, i.e., 
array, the "id" information of each item in sequence list is implied, and which 
will be used later. Just now, we discussed the details of two sequence lists, in 
the sequel, we present the details of two labels. 

Two labels: As we know, in the standard plane sweep method, when a point p 
(G Q) is a left endpoint of certain segment, we need insert the segment into the 
red-black tree R. In our proposed method, we assign two labels, say lb\ and Ibz, 
to the segment Q before it is inserted into the red-black tree R. 

The first "label" is Boolean type, which identifies that a segment is coming 
from which one circular-arc polygon. Specifically, the value of lb\ is "true" if the 
segment from CAP%, otherwise, the value of lb\ is "false". Recall standard plane 
sweep method, when two segments are adjacent, then checking if they intersect 
or not. In our proposed method, if the first labels of two adjacent edges have 
same value, we do not check if they intersect or not. 

For example, there are a number four of processed related edges in Figure 
[31 we assign "true" to the first labels of ab and 6c, and assign "false" to the 
first labels of AB and BC. At the second sweep line (see Figure [3(b) [ ), edge be is 
inserted into red-black tree R, in here, though edge ab immediately below the new 
inserted edge be, we do not check if they intersect or not, since they have same 

3 

Note that, the segment discussed in here is refer to the processed related edge 
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value in terms of the first labels. By doing so, it not only can avoid checking 
if they intersect each other, but also can avoid the vertexes from circular-arc 
polygons being reported as intersection points. 

The second "label" is an integer type denoting a serial number, which corre- 
sponds to the "id" of item in sequence list. Regarding to the example in Figure 
[3J we assign "1" to the second label of edge ab, and assign "2" to the second 
label of edge be. Similarly, we assign "1" and "2" to the second labels of edge 
AB and BC, respectively. By doing so, once we detected an intersection point, 
we, based on their labels, can quickly find the item stored in sequence list. Then, 
to insert a pair of decimals to denote the locations of intersection points. 

For example, see Figure [3(b)| when we detected the intersection point I\, 
since the lb% and Ibi of edge ab are "true" and "1" , respectively, therefore, we 
can quickly know that we should insert a decimal into the iteml of sequence list 
Si. Similarly, since the lb\ and Ib-i of edge BC are "false" and "2", respectively, 
we also can quickly know that we should insert a decimal into the item2 of 
sequence list Si- 

Otherwise (i.e., if we do not adopt the two labels), once we detected an 
intersection point, we have to scan the sequence list in order to insert the decimal 
in appropriate item, this manner is inefficient especially when the number of 
items in the two sequence lists are much more. 

It is worth noting that the tactics discussed just now also are effective for 
dealing with Boolean operation on traditional polygons. Up to now, we have 
discussed the main tactics used for constructing two sequence lists. In the next 
subsection, we present the specific algorithm. 

Algorithm for constructing two sequence lists Figure [B] illustrates the 
algorithm for constructing the two sequence lists. The general idea is to choose 
those related edges based on the extended boundary lines at first (line 1-2). 
Next, We construct two empty sequence lists and initialize them (line 3-5). In 
the process of initializing the two sequence lists, non-x-monotone arcs are de- 
composed, then processed related edges are put into the sequence lists in order, 
and appropriate values are assigned to the tri- value switchs (line 25-34). After 
this, we compute the intersection points (line 6-23). Note that, in the process 
of obtaining the intersection points, two "labels" are assigned to the segment 
before it is inserted into red- black tree (line 10), and we use the two sequence 
lists to store the locations of intersection points (line 13, 18, 23), by doing so, it 
can let these intersection points be well organized. 

Theorem 2. Given two circular- arc polygons with M and N edges, respectively, 
and assume there are I related edges between the two polygons, then, to construct 
the two sequence lists takes 0(M + N + I + (I + k) * logZ) time, where k is the 
number of intersections. 

Proof Choosing related edges needs to find the MBRs of the two circular- 
arc polygons, and to determine the effective axis, by comparing the horizontal 
and vertical spans of the intersection set of two MBRs. Next, based on two 
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Procedure Constructing_Two_Sequence_l_ists 
Input: CAPx and CAP 2 

Output: Si and S2, Ri and R2 // the locations of intersections has been 

/ / well organized in Si and S2 

(1) Find the MBRs; to determine the effective axis and the extended boundary lines 

(2) R\ <— related edges from CAP±; R2 <— related edges from CAP2 // Dcfinition[6l 

(3) Create two empty sequence list Si and S 2 

(4) for each i £ [1, 2] 

(5) Subprocedure lnitilize_Sequence_List (Ri, Si) 

(6) Sort the endpoints of all the segments from Si and S2 (i-e., processed related edges), and 
store all the endpoints in a heap Q 

(7) H 4— / / R is a red-black tree used for maintaining the order of segments 

(8) for each point p in Q 

(9) if (p is a left endpoint of segment s) 

(10) Assign two "labels" to s, and insert s into R. 

(11) if ( s and the segment immediately below or above it are not from same polygon ) 

(12) if (s intersects the segment) 

(13) Insert the intersection point into Q and insert a pair of decimals into Si and S2; 

(14) else if (p is a right endpoint of segment s) 

(15) if ( the pair of segments directly above and below s are not from same polygon ) 

(16) if (there exist intersection point between them and this intersection point Q ) 

(17) Insert the intersection point into Q; delete s from R 

(18) Insert a pair of decimals into Si and S2; 

(19) else // p is an intersection point of two segments, say s and t 

(20) Swap the position of s and t in R; / / assume s is the upper segment now 

(21) if (s/t and the segment above/below it are not from same polygon) 

(22) if (they intersect each other) 

(23) Insert the intersection point in Q and insert a pair of decimals into Si and S2; 

(24) return Si and S2 , Ri and R 2 
Subprocedure lnitilize_Sequence_List 
Input: Ri, Si 

(25) temp <— 1 // temp is used for set the tri-value switch 

(26) for each related edge £ Ri 

(27) if (it is a non-x-monotone circular arc ) 

(28) Decompose it and put the decomposed arcs into Si, each tri-value switch 
is set to "iemp" 

(29) if (temp=l) 

(30) temp ^2; 

(31) else / / temp— 2 

(32) temp 

(33) else // it is not a non-x-monotone circular arc 

(34) put it into Si, in which the tri-value switch is set to "0" 

Fig. 6. The Algorithm of Constructing Two Sequence Lists 



extended boundary lines, we compare the geometry relation between each edge 
and extended boundary lines. So, line 1-2 takes 0(M + N) time. 

Creating two empty sequence lists takes constant time. In order to initialize 
the two sequence list, we need to decompose each non-x-monotone arc, which 
takes constant time, at worst case, all the related edges are non-x-monotone arcs. 
In addition, according to Theorem [TJ there are no more than 3 * / items in the 
two sequence lists, so initializing two sequence lists takes 0(1) time. Sorting all 
endpoints of segments (from Si and S2) in heap Q takes 0(1 * log/) time. In 
addition, initializing red-black tree R takes constant time. Thus, line 3-7 takes 
0(Z + J*logZ) time. 

Since there is are no more than 3 * / segments in Si and 52, the number of 
endpoints of segments are no more than 6 * /, thus, the number of times the 
for loop (line 8-23) is executed is no more than 6 * I + k. Since the number of 
segments in R never exceed 3 * Z, each operation (e.g., insert, delete, swap, find 
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above/below segment) on R within the for loop can be performed in O(logZ) 
time. In addition, all the other operations within the for loop can be performed 
in constant time, these operations include: (1) assigning segment s with two 
labels; (2) determining if two segments are from same polygon; (3) determining 
if two segments intersect each other; (4) inserting a pair of decimals into Si and 
6*2; (5) determining if a point p (in Q) is a left/right endpoint of segment s, 
or an intersection point. Thus, line 8-23 takes 0((6 * I + k) * log/) time, i.e., 
0((l + k) * log/) time. 

At last, returning S\ and S2 takes constant time. Putting all together, we have 
that the time for constructing the two sequence lists is 0(M+N+l + (l+k)*\ogl). 
□ 

3.2 Constructing the modified polygons 

In previous subsection, we discuss how to obtain two sequence lists, in which the 
locations of intersection points are well organized. In this subsection, we discuss 
how to build two modified polygons (i.e., two new linked lists) that are used for 
traversing. 

Since the two sequence lists obtained in previous step are built along the 
direction of counter-clockwise w.r.t. original circular-arc polygons, this makes 
the process of building two new linked lists rather simple. What we need to do is 
using the information stored in two sequence lists to replace those related edges 
in original circular-arc polygons. Recall that the original circular-arc polygons 
also are double linked lists, we have discussed the data structure in subsection 

i2~n 

Note that, there are three important issues that we have to take them into 
consideration. 

1. In the previous step, we use decimals to denote the locations of intersec- 
tion points, therefore, we have to transform a series of decimals into the 
coordinates of intersection points. 

2. If there exist intersection points on arcs, we have to insert new appendix 
points into these arcs, which is ready for eliminating ambiguity. 

3. We decomposed non-x-monotone arcs into x-monotone arcs ever, then, we 
should combine them in the process of replacing. 

In the sequel, we first discuss the details of the above three issues, followed 
by presenting the specific algorithm for constructing two new linked lists. 

Transforming the decimal into the coordinate Given any item in sequence 
list Si or S2, if the second domain is not null, it indicates that there is at least one 
intersection point on this segment. In the process of replacing, we have to trans- 
form decimals into the coordinates of intersections, the process of transforming 
from a decimal into the coordinate of intersection point is straightforward, it can 
be summarized as follows. 
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case 1: If the item is a straight line segment, we have 



[xi = (x 2 -x 1 )*S + Xi ^ ^ 

\yi = (y2-yi)*$ + yi 

where 8 is the decimal, (xi,yi) is the coordinate of intersection point, (xi,yi) 
and (£2,2/2) are the coordinates of two endpoints of the segment, respectively. 

case 2: If the item is a circular arc, say abc. We first compute the centre and 
radium, denoted by {xo,yo) and r respectively, of circle corresponding to the 
circular arc abc, as shown in Figure [7] Next, we compute the central angle, 
denoted by 9, of aJc based cosine theorem. Since & = arc \ engt t h h f f rom a \° 1 = 6, 

J 1 6 arc length from a to c 1 

thus, we have 9\ =9*6. Let us consider polar plane, since the coordinate of 
endpoint a is given, further, since x = xq + r * cosa and y = yo + r * sina, thus, 
we can get 02- At last, we have 

\xi = xo + r*cos(9*S + 9 2 ) 
\yi =yo + r * sin(9 * 5 + 9 2 ) 

where 6 is the decimal, (xi,yi) is the coordinate of intersection point. 




Fig. 7. Example of Transforming from Decimal to Coordinate 



Inserting new appendix points In subsection 12.11 we introduce the data 
structure for representing the circular-arc polygon, in which we always add an 
appendix point between two vertexes if the edge is a circular arc rather than a line 
segment. In the process of replacing, we also have to ensure this property if the 
edge is a circular arc. Note that, when the intersection points appear on a circular 
arc, this arc will be decomposed by these intersection points, therefore, for each 
sub-segment, we have to ensure there is an appendix points for eliminating the 
ambiguity. So, we have 

Theorem 3. Given there are a number k of intersections on a circular arc, then 
we need insert at least k and at most k + 1 new appendix points for eliminating 
the ambiguity. 
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Proof Since a number k of intersections can subdivide a complete circular arc 
into k + 1 small circular arcs, and for each small circular arc, one appendix point 
is needed and is enough to eliminate the ambiguity. So, for k + 1 small circular 
arcs, k + 1 appendix points are needed. Further, since there is an appendix point 
beforehand. Therefore, when there is no intersection point coinciding with this 
appendix point, only k new appendix points are needed. Otherwise, k + 1 new 
appendix points are needed. □ 

Merging the decomposed arcs In order to merge the decomposed arcs, a 
bruteforcing solution is that, after we obtained the resultant polygon, we then 
check each pair of adjacent edges, to see if they can be merged. Obviously, this 
way is inefficient. 

In here, we utilize the information stored in tri-value switch to help us up- 
grade the efficiency of merging the decomposed arcs. Specifically, when we de- 
tected that the value of tri-value switch of an item is "1" or "2" , we continue to 
fetch its next item from sequence list, in this way, until a group of consecutive 
items, which have same value in terms of the tri-value switch, are fetched from 
sequence list. Then, we do following processing. 

— If the second domain of each item is null, we discard the fetched items. This 
is because there is no intersection point on these decomposed arcs, the result 
of merging them should be as same as the original edge. 

— Otherwise, we transform the decimals into the coordinates of intersection 
points based on the method discussed earlier, and insert new appendix point 
for eliminating ambiguity, then merge them and replace the edge in original 
circular-arc polygon. 

As an example, let us revisit Figure [SJ there are a number eight of items in 
sequence list Si, the values in their tri-value switches are "0", "1", "1", "2", "2", 
"1" , "1" , "0" , respectively. Assume the first item has been processed, then, for 
the second item, the value of its tri-value switch is "1" , in this case, we continue 
to fetch its next item, i.e., the third item. After this, for the fourth item, the 
value of tri-value switch is not "1" , therefore, we stop fetching this item for 
the present. Note that, both of the second domains of the second and third 
item are null, so we discard the two items rather than combing them. Similarly, 
for the fourth and fifth items, we also discard them. For the sixth and seventh 
items, they have same value in terms of their tri-value switches, and there is a 
decimal in the seventh item, thus, we need transform the decimal into coordinate 
of intersection point, and insert new appendix point for eliminating ambiguity. 
More important, we need merge the two items (i.e., two decomposed arcs), and 
use the merging result to replace the edge in original circular-arc polygon. 

Algorithm for constructing two new linked lists Figure [5] depicts the 
algorithm for constructing two new linked list. The general idea is to judge, for 
each edge in original circular polygon, if it is a related edge. If it is not (line 
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Procedure Const ructing_Two_New_Linked_Lists 

Input: CAPx and CAP 2 , S\ and S 2 , Ri and R 2 
Output: CAP? and CAP 2 * 

(1) Initialize two empty linked list CAP? and CAP£ 

(2) for each % £ [1,2] 



(3) id <— 1 1 1 id corresponding to the serial number of item in S.% 

(4) for each edge e in CAPi / / Recall our data structure, two or three nodes denote an edge 

(5) if (e g Ri ) 

(6) Add e to CAP? 

(7) else // e is not a related edge 

(8) tri the value of tri-value switch in Si[id] / j Si[id] denotes the id 1 item in Si 

(9) if (tri— 0) // not a decomposed arc 

(10) if (the second domain of Si[id] is null) //there is no intersection point 

(11) id <~ id+ 1; add e to CAP* 

(12) else // the second domain of Si[id] is not null 

(13) if (Si[id] is a circular arc) 

(14) Transform the decimal into coordinate // Equation |3.2| 

(15) Insert new appendix point for eliminating ambiguity 

(16) else // Si[id] is a straight line segment 

(17) Transform the decimal into coordinate // Equation l3.il 

(18) Add the information transformed from Si[id] into CAP?; id <— id + 1 

(19) else / / tri is equal to 1 or 2. 

(20) j <- 

(21) repeat / / copy the consecutive decomposed arcs 

(22) j <- j + 1; tempp] <- Si [id] 

(23) id 4— id + 1; tri* -f— the value of tri-value switch in Si [id] 

(24) until tri ^ tri" 

(25) if (the second domain of each temp[j] is null) // except j — 

(26) Add e to CAP? 

(27) else 

(28) mer Result <— //mer Result is used for storing the merging result 

(29) for each j // except j — 

(30) if (the second domain of ternp[j] is not null) 

(31) Transform the decimal into coordinate of intersection / /Equation l3.2"1 

(32) Insert new appendix point for eliminating ambiguity 

(33) raerResult <— Merge these transformed information 

(34) Add combResult into CAP" 

(35) return CAP? and CAP 2 * 



Fig. 8. The Algorithm for Constructing Two New Linked Lists 



7-30), we further differentiate if it is a decomposed arc, this is easy to confirm 
based on the value of tri-value switch (line 8). Line 9-18 is used for dealing with 
the case where the edge e is not a decomposed arc. To contrary, line 19-34 is 
used for handling the case where the edge e is a decomposed arc. When it is a 
decomposed arc, we first fetch all the consecutive decomposed arcs (line 20-24), 
then to see if there is intersection point on these decomposed arcs. If it is not, we 
need not merge them (25-26). Otherwise, we merge them and add the merging 
result into new linked list (line 27-34). 

Theorem 4. Given we have obtained the two sequence lists Si and S2, then, to 
construct the two new double linked lists (CAP^ and CAP^ ) takes 0(M + N + 
I + k) time. 

Proof Transforming each decimal into coordinate of intersection point takes 
constant time, there are k intersections, so it takes 0{k) time for transforming. 
Inserting each appendix point takes constant time, even if all the intersection 
points locate on arcs rather than on line segments, according to Theorem [31 
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there are no more than 2 * k new appendix points, so inserting appendix points 
also takes 0(k) time. 

Since the number of edges in CAP\ and CAP 2 is M + N, the number of 
times of the second for loop is executed is M + N. More specifically, the number 
of times of Line 6 is executed is (M + N — I), the number of times of line 7- 
34 is executed is I. Note that, we no longer consider the time for transforming 
decimals to coordinates, and the time for inserting new appendix points, since 
we have analysed them in previous paragraph. At worst case, each related edge 
is non-x- monotone arc, according to Theorem [TJ so the number of times of line 
21-24 is executed is no more than 3 * I, further, since each single operation in 
line 7-34 takes constant time, therefore, line 7-34 takes 0(3 * I) time. 

Thus, we have that constructing the two new double linked lists takes 0(k + 
k + M + N-l + 3*l) time, i.e., 0(M + N + k + I) time, when Si and S 2 is 
given beforehand. □ 

3.3 Assigning entry-exit property and traversing 

In this subsection, we first discuss how to assign entry or exit property to inter- 
section points, followed by presenting three traversing rules, which are used for 
obtaining the intersection, union and difference, respectively. 

Entry-exit property Once we obtain two new linked lists, then, we assign the 
intersections with entry or exit property alternately. For ease of discussion, we 
use a simple example, in which there are less edges. 




(a) 



iteml 

itme2 
item3 
itme4 



1 (V1V2) 


0.25,0.7 fo] 




(V,V 3 ) 


NULL | 0~| 



(V 3 V 4 V 5 ) 


0.2 fo] 


1 (V5V0 


0.15 fo] 



iteml 
itme2 
item 3 
itme4 



(V 6 Vy) I I NULL I [0] 



(V 7 V S ) I 1 0-2.0-5 I \T\ 



(V s V,Vio)| [null] \T\ 



(ViqV<,) I I 0.2,0.6 I \J} 



(b) 



Fig. 9. Example of Two Simple Circular-arc Polygons and Their Sequence Lists 



Taking the two circular-arc polygons in Figure 9(a)| as an example, CAP\ = 
{VI, V 2 , V3, V4, V5} and CAP 2 — { Vg , Vj,Vs,Vt), Vw}, we assume we have ob- 
tained the two sequence lists Si and S2 (see Figure |9(b)| based on the al- 
gorithm in Figure |6] Further, it is easy to obtain the two new linked lists 
CAP* and CAP^ using the algorithm in Figure [H More specifically, CAP{ = 

{Vi,I 2 ,h, V 2 , V 3 ,Ni,h,% V 5 ,h}, and CAP 2 * = {V 6 , V r , J 3 , h, V & ,V 9 , V 10 , h,I 2 }. 
Then, for the new linked list CAP^, we assign ll entry, exit, entry, exit" to il I 2 , 
I3, I4, Ji", respectively. The entry-exit properties of " I 2 , I3, I4, /1" in CAP£ 
are as same as that in CAP{ . The results are illustrated in Figure [TUJ 
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head 



entry exit 



entry 



exit 



V, <-> h *-*■ h <-* V 2 ■*->• V 3 <-» Ni "*-*■ I 4 V 4 <-> V 5 - * h 



head 



(a) CAP* 
exit entry 



T 



v 7 



I 3 



I4 



V, 



V, 



exit entry 



Ii 



I2 



(b) cap; 

Fig. 10. Example of Entry and Exit Properties 



Traversing the two new linked lists Once the entry-exist properties are 
assigned to intersection points, we then traverse the two new linked lists for 
obtaining the resultant polygon. There are some differences among three opera- 
tions (i.e., intersection, union and difference) when we execute traversing. In the 
sequel, we discuss the three operations one by one. 

Intersection: We first choose an intersection point with entry property from 
CAPl as the starting point, we then traverse along the direction of counter- 
clockwise, once we meet an intersection point with exit property, we shift to 
CAP^, and traverse along the direction of counter-clockwise. Similarly, if we 
meet an intersection point with entry property in CAP 2 *, then shift back to 
CAP*. In this way, a circuit, say CC%, will be produced. After this, we check if 
there is another intersection point with entry property that is not a vertex of the 
produced circuit. If no such intersection point, then we terminate traversing, the 
CC\ will be the intersection between the two original circular-arc polygons CAP\ 
and CAP2. Otherwise, if there is still intersection point with entry property that 
is not a vertex of any produced circuit, we let the intersection point as a new 
starting point, and to traverse using same method discussed just now, until no 
such intersect point exists. In the end, we get multiple circuits, which are the 
intersection between the two circular-arc polygons. 




Fig. 11. Example of Intersection 



To continue previous example, we first choose I2 in CAP* as the starting 
point, as shown in Figure [TT1 Next, since I 3 has exit property, we shift to CAP^. 
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Similarly, since J4 has entry property, we shift back to CAPf. In this way, a 
circuit CC\ is produced (see the dotted lines). Moreover, there is no other inter- 
section point^with entry property that is not in CC\. Therefore, the intersection 
is {I 2 ,h,h,V 4 ,V 5 ,h}. 

Union: We first choose a vertex from CAPf such that it does not locate in 
CAP 2 , as the starting point. Then, we traverse along the direction of counter- 
clockwise, once we meet an intersection point with entry property, then we shift 
to CAP 2 *, and traverse along the direction of counter-clockwise. Similarly, if we 
meet an intersection point with exit property in CAP 2 , then shift back to CAP£ . 
In this way, a circuit will be produced, this is just the union between the two 
original circular-arc polygons CAP\ and CAP 2 . 




Fig. 12. Example of Union 



Regarding to previous example, in order to get the union, we first choose 
Vi as the starting point, as shown in Figure [12] Next, since I 2 has entry prop- 
erty, we shift to CAP 2 . Similarly, since I3 has exit property, we shift back to 
CAP*. In this way, until we are back to the starting point V\. The circuit 
(see the dotted lines) is just the union between CAP\ and CAP2, which is 
{Vx,I 2 ,V 6 , V 7 , h,V 2 ,V 3 , N u h, V 8 ,%, V w , h}. 

Difference: We first choose a vertex from CAPf such that it does not lo- 
cate in CAP 2 , as the starting point. Then, we traverse along the direction of 
counter-clockwise , once we meet an intersection point with entry property, then 
we shift to CAP 2 *, and traverse along the direction of clockwise . Similarly, if 
we meet an intersection point with exit property in CAP 2 , then shift back to 
CAPi. In this way, a circuit, say CC\, will be produced. In addition, we check 
if there is another vertex such that (1) it is a vertex of CAPf; (2) it is not a 
vertex of any produced circuit; (3) it does not locate in CAP 2 * . If no such ver- 
tex, then we terminate traversing, the GC\ will be the difference between the 
two circular-arc polygons. Otherwise, if there is still vertex that satisfies the 3 
conditions above, we let the vertex as a new starting point, and to traverse using 
same method discussed just now, until no such vertex exists. In the end, we get 
multiple circuits, which are the difference between the two circular-arc polygons. 

For instance, in order to get the difference between previous two circular- 
arc polygons, we first choose V\ as the starting point, and begin to traverse 
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along the direction of counter-clockwise. Since I2 has entry property, we shift 
to CAP2, note that, we traverse along the direction of clockwise in CAP 2 *, 
therefore, the next node of 1% is I\. Since I\ has exit property, we shift back to 
CAPi . In this way, a circuit CC% (see black dotted lines) is produced, as shown 
in Figure[l3l After this, we choose V2 as a new starting point since it satisfies the 
three conditions above. Then, we begin to traverse using same method discussed 
just now, we get another circuit CC2 (see the grey dotted lines). Thus, the 
difference between CAP\ and CAP 2 consists of two parts, i.e., {Vi,l2, h} and 

Note that, though the traversing rules have some minor differences among 
three operations, their time complexity are equal, we in here only analyse the 
time complexity on intersection operation for saving space, regarding to other 
two operations, they can be derived similarly. 

Theorem 5. Given the two new double linked lists CAP* and CAP 2 , to obtain 
the resultant polygon takes 0(k + m + n + I) time. 

Proof Assigning entry-exit property to each intersection point takes constant 
time, there are k intersections on each new double linked list, so it takes O(k) 
time. 

There are no more than "M + N + 3 * /" edges in CAP* and CAP 2 , so the 
number of times for visiting each node in the two new double linked lists is no 
more than M + N + 3 * I, each operation on a node (e.g., determining the type of 
a node, inserting a node to resultant polygon) takes constant time, so traversing 
takes 0(M + N + I) time. 

Pulling the above two results together, thus, we have that to obtain the 
resultant polygon takes 0(M + N + k + I) time, when CAP{ and CAP 2 are 
given beforehand. □ 



3.4 The total running time 

Up to now, we have addressed all steps for achieving Boolean operation on 
circular-arc polygons, and analysed the running time of each step. In this sub- 
section, we give the total time complexity of our propose approach. 



22 



The total time complexity can be derived by pulling the results in Theorem 
GO H and [5] together. 

CBoolean =0( (M + N + I + {I + k) * \ g I) + (M + N + k + I) + (M + N + k + I) ) 

=0{ M + N + (I + k) * (1 + logi) ) 
=0( M + N+ (l + k) *logZ ) 

(3.3) 

where C Boolean is the total complexity. Thus, we have the following Theorem. 

Theorem 6. Given two circular-arc polygons with M and N edges, respectively, 
and assume there are a number I of related edges between the two circular- arc 
polygons. Then, to achieve the Boolean operation between them takes 0(M + 
N + (/ + k) * log?) time, where k is the number of intersections. 

4 Comparison 

Though the theme of this paper is Boolean operation on circular-arc polygons, 
our proposed method can be easily extended to deal with Boolean operation on 
traditional polygons. In particular, it is easy to verify that the running time for 
Boolean operation on traditional polygon also is no more than 0(M+N+(l+k)* 
logZ), this improves state-of-the-art results, since "/ < M + N" is always true. 
Moreover, as stated in Section [TJ though Gong et al. [19] proposed algorithm 
also can deal with boolean operation on circular-arc polygons, but the execution 
time required by their algorithm is 0(M * N), which is equal to the execution 
time required by the naive method 1. 

For ease of comparison, we summarize our results (including naive methods) 
and several state-of-the-art results, which are related to this work. Note that, 
in some papers, the authors adopted 'n' to denote the number of edges of two 
polygons. For clarity, we rewrite their results in form of M and N, where M+N = 
n. These results are listed in Tabled] 

Table 1. Results comparison 



Method 


Description 


Time complexity 


Martinez et al. [5] 


Traditional polygon, 2009 


0{ (M + N + fc) * log(M + N) )=0(M + N + k + (M + 






N + k) * log(M + TV)) 


Liu et al. [9] 


Traditional polygon, 2007 


0(M + N + k + (M + N + k) * log(Af + N)) 


Peng et al. fill 


Traditional polygon, 2005 


0{ M + N + k + (M + N) * (k + log(M + TV)) + M * N) 


Rivero et al. 1 101 


Traditional polygon, 2000 


0( (M + k) * (N + k) ) 


Grciner et al. \7\ 


Traditional polygon, 1998 


0( M * N ) 


Gong ct al. EH 


Conic polygon, 2009 


0(M * AT) 


Naive approach 1 


Circular-arc polygon, - 


0( M * N ) 


Naive approach 2 


Circular-arc polygon, - 


0( (M + N + k) * (M + N + log(M + N) ) ) 


The proposed method Circular-arc polygon, - 


0(M + N + (/ + k) * log/) 
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5 Conclusion 



This paper investigates Boolean operation on circular-arc polygons. By well con- 
sidering the nature of problem, efficient data structure and targeted algorithms 
are proposed. The data structure adopts a double linked list, in which the ap- 
pendix point is designed for eliminating the ambiguity. In the proposed algorithm, 
several tactics are presented for enhancing the efficiency, the running time of the 
proposed algorithm is 0(M + N + (I + k) * log/). On the whole, the proposed 
method is concise, easy to understand and implement, but without loss of the 
efficiency. 

Note that, the following situations are possible when we consider Boolean 
operation on circular-arc polygons: (1) there possible are some degenerations, 
e.g., vertex of circular-arc polygon CAP\ touches the edge of circular-arc polygon 
CAPi\ and (2) there possible are many holes in circular-arc polygon; and (3) the 
circular-arc polygon possible is self-intersection. Throughout this paper, these 
problems are not addressed, we leave them as our future work. 
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